Critical Scaling at Zero Virtuality in QCD 
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We show that at the critical point of chiral random ma- 
trix models, novel scaling laws for the inverse moments of the 
eigenvalues are expected. We evaluate explicitly the pertinent 
microscopic spectral density, and find it in agreement with nu- 
merical calculations. We suggest that similar sum rules are of 
relevance to QCD at the critical temperature, and even above 
if the transition is amenable to a Ginzburg-Landau descrip- 
tion. 

PACS numbers : 11.30.Rd, 11.38.Aw, 64.60. Fr 



1. A large number of physical phenomena can be 
modeled using random matrix models [0. An impor- 
tant aspect of these models is their ability to capture the 
generic form of spectral correlations in the ergodic regime 
of quantum systems. This regime is reached by electrons 
traveling a long time in disordered metallic grains Q or 
virtual quarks moving a long proper time in a small Eu- 
clidean volume |||. 

In QCD, the ergodic regime is characterized by a huge 
accumulation of quarks eigenvalues near zero virtuality. 
This is best captured by the Banks-Casher Q relation 
\{qq) \ = £ = irp(0), where the nonvanishing of the chiral 
condensate in the vacuum signals a finite quark density 
p(X = 0) at zero virtuality. This behavior is at 
the origin of spectral sum rules J5| , which are reproduced 
by chiral random matrix models ||. These sum rules 
reflect on the distribution of quark eigenvalues and cor- 
relations 0. 

If QCD is to undergo a second or higher order chiral 
transition, then at the critical point there is a dramatic 
reorganization of the light quark states near zero virtu- 
ality as the quark condensate vanishes. In section 2, we 
suggest that such a reorganization is followed by new scal- 
ing laws, which are captured by a novel microscopic limit. 
In section 3, we use a chiral random matrix model with 
a mean- field transition to illustrate our point. In section 
4, we explicitly construct the pertinent microscopic spec- 
tral distribution in the quenched case and compare it to 
numerical calculations. In sections 5,6 we argue that the 
sum rules following from the new scaling, may be of rel- 
evance to QCD at the chiral critical point, provided that 
the transition follows the general lore of universality, and 
suggest that spectral correlations persist even above the 
transition temperature. 

2. The nonvanishing of £ in the QCD vacuum im- 
plies that the number of quark states in a volume V, 



N(E) = V J dXp(\) in the virtuality band E around 
grows linearly with E, that is N(E) ~ EV. As a re- 
sult, the level spacing A = dE/dN ~ 1/V for N ~ 1, 
and the eigenvalues of the quark operator obey spectral 
sum rules (^]. During a second or higher order phase 
transition £ vanishes in the chiral limit. Scaling argu- 
ments give £ ~ m 1 / 5 at the critical point, where m 
is the current quark mass 0. It follows again from 
the Banks-Casher relation [[l|, that for small virtualities 
A, ~ IAI 1 / 5 to leading order in the current quark 

mass §. Hence, N(E) ~ VE 1+1 / S , and the level spac- 
ing is now A* ~ l/V s ^ s+1 '> at iV ~ 1. 

For a mean-field exponent S = 3, and we have A* = 
V~ 3 / 4 , which is intermediate between V~ x in the spon- 
taneously broken phase and T/ -1 / 4 in free space. At 
the critical point there are still level correlations in the 
quark spectrum except for the free limit, corresponding 
formally to S = 1/3. We now conjecture that at the crit- 
ical point, the rescaling of the quark eigenvalues through 
A — * A/ A*, yields new spectral sum rules much like the 
rescaling with A = 1/V in the vacuum Q. The mas- 
ter formula for the diagonal sum rules is given by the 
(dimensionless) microscopic density of states 



i/*(s) = lim (VA* 



p(sA*) 



(1) 



and similarly for the off-diagonal sum rules in terms of 
the microscopic multi-level correlators. Since we lack an 
accurate effective action formulation of QCD at T = T c 
(a possibility based on universality is discussed below), 
the nature and character of these sum rules is not a pri- 
ori known, but could easily be established using lattice 
simulations in QCD. 

Could these sum rules be shared by random matrix 
models? We will postpone the answer to this question 
till section 5, and instead show in what follows that the 
present scaling laws hold at zero virtuality for chiral ran- 
dom matrix models with mean-field exponents. 



3. Consider the set of chiral random plus deterministic 
matrices 

im t + A 
t + A^ im 



M = 



(2) 



where 4 is an W x Af complex matrix with Gaussian 
weight, m a 'mass' parameter, and t a 'temperature' pa- 
rameter. Such matrices or variant thereof have been in- 
vestigated by a number of authors in the recent past [|l0| . 
Their associated density of states is 



1 



p(X) = —{TrS(\-M)) 



(3) 



where the averaging is carried using the weight associated 
to the following partition function 



Z[m,t] = J dAdA^ detM N ?e- NTrAAf . 



(4) 



For t = m = the density of states is p(Q) = 1/V, while 
zero for t > 1 and m = 0. At t = 1, /d(A) = IAI 1 / 3 , which 
indicates that the i-driven transition is mean-field with 
p(0) as an order parameter. In particular, at t = 1 the 
level spacing is A* = N~ 3 / 4 near zero. 

Standard bosonization of the partition function (|4|) 
yields 

Z[m, t] = JdPdP^ e ^o S det[ [m+ p)(rM-pi) + t°]-NTrppi (5) 

where P is an Nt x iV/ complex matrix. We may shift P 
by the mass matrix P — > Q = P + m and get 



Z[m, t] = J dQdQ^ det(QQ f + i 2 )^ 

Xe -iv[TrQQt-mTr(Q+Q t )+m 2 ] 



(6) 



dropping the irrelevant normalization factor. We will 
now specialize to the critical temperature t = 1, and 
denote the rescaled mass by x — im/A, and eigenvalue 
by s = A/ A*. This suggests the rescaling Q — > Q = 
AT 1 /4q ! go that 



= /dQdQte-iTrfQQtj^^TrW+Qt 



reducing to 



Z[x) 



rdr e 



(7) 



(8) 



n 



for one flavor. Expanding this integral in powers of x, 

l k (m/A*) 2k f k + 1 
,, =1 (fc!) 2 ^ V 2 



Z[m] = Z[0] 



(9) 



gives rise to spectral sum rules for the moments of the 
reciprocals of the eigenvalues of M by matching the mass 
power in the spectral representation of the partition func- 
tion, 



Z[m]/Z[0] = (l[ fl + ^T) 

\A fc >0 V A kJ 



(10) 



where the averaging is done over the Gaussian random- 
ness with the additional measure nA fc >o^fe- Matching 
the terms of order m 2 yields 



y — 



\A fc >0 



(11) 



and matching the terms of order m 4 gives 



^ A? 

\A fc >0 k , 



At 



\A fc >0 



The relations ( |ll[]l2| ) are examples of microscopic sum 
rules at the critical point t = 1. 

One should note here that the preceding calculation 
has been performed for the gaussian matrix model. It 
turns out that if we were to add terms of the form 



-JV ff oTr(PP + ) 2 -Ar £ ,i(TrPP t 



(13) 



to the measure in (||), then the spectral sum rules are 
in general unchanged. Indeed, for Nf = 1 the addition 
amounts to a global shift x — » x/(l + 2go + 2gi). Higher 
powers of PP* are subleading after rescaling and do not 
affect the sum rules in large N. 

4. The diagonal moments of the reciprocals of the 
rescaled eigenvalues are generated by the microscopic 
density (Q), in the limit N — > oo and A — * but s = A/ A* 
fixed. The microscopic density (Q) for Nf flavours in a 
fixed topological sector n : v*,N f .n, can be evaluated us- 
ing supersymmetric methods [p]Jll[. For example, 



f*,o,o(s) 
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.Ms) Jo (s) + k Q (s)j 2 (s) - fci (s)j'i (s)) 



(14) 



where k n (s) and j n (s) are given by 



Us) = ^ dz Z 2n+lj o(2z _±_) e -^/2 (15) 



k n (s) = [dz z ^+2 Kl( 2"y 4 



/2 



(16) 



and where the integration contour C is the sum of two 
lines: C = [— (1 + i)oo, 0] U [0, (1 — i)oo]. In this spirit, 
the first sum rule (0) reads 



1 f u 2 
-j- ^*,i,o(s) ds = —= 



(17) 



The second sum rule (^2|) involves the 2-level microscopic 
correlator for Nf = 1 and n = which can be obtained 
using a similar reasoning. 

In Fig. |l| we compare the quenched Nf — 1 (left) and 
unquenched Nf = 1 (right) results to the numerically 
generated microscopic spectral density at t = 1 using 
N = 100 size matrices. The agreement suggests that 
the present method of finding the scaling properties of 
microscopic spectral distributions can be used to accu- 
rately determine the value of the critical exponent 6 in 
lattice simulations. 
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FIG. 1. left: Vt,o,o(s) at t = 1 and Nf = for matrices of 
size N=100 (dots) and the theoretical prediction (^) (solid 
line), right: f*,i,o(s) at t = 1 and iV/ = 1 for matrices of size 
N = 20 (dotted), iV = 50 (dashed) and JV = 100 (solid). 



5. In QCD the finite temperature transition is still 
debated. For very light quark masses, it is suspected 
to be second order, although a cross over is not ruled 
out. The character of the transition depends crucially 
on the number of flavors Nf and the fate of the ?7a(1) 
quantum breaking. For two light flavors, we may assume 
with Pisarski and Wilczek [jl2j that the transition is from 
an SU(2) spontaneously broken phase to Z2 x SU(2) x 
SU(2) ~ 0(4) |l3[ ]. Choosing the order parameter <f> = 
a + it ■ 7r (vacuum analogous to a ferromagnet), implies 
for the Ginzburg-Landau potential 



V($) = +m Tr($ t + $) 
+. 9o (T)Tr($t$) 



5l (T)(Tr($t$))2 



(18) 



at zero vacuum angle 9. The dots refer to marginal or 
irrelevant terms. For a second order transition, go{T) ~ 
T — T c , which is negative below T c and positive above. 
Note, that in this section we analyze this hamiltonian 
just from the mean-field point of view by truncating it to 
the space of constant modes. 

For T = T C the potential |[|) is reminiscent of the one 
for the chiral random matrix model discussed above 
If we note that the measure on the manifold with re- 



stored symmetry is 



-I3V 3 V 



-vv 



we conclude that 



( |18| ) is enough to accommodate for the level spacing 
A* = 1/V 5 I (5+ V with 6 = 3 (mean-field). Indeed after 
the rescaling x = im/A* we recover (Q) with g\ = 1/2 
and the proper identification of the manifold. 

For T > T c we can use (|l^) to define new sum rules 
for the quark eigenvalues in the sector with zero winding 
number ||. In particular (Nf — 2) 




(19) 



The rhs measures the variance in the scalar direction on 
an invariant 0(4) manifold with e~ vv as a measure. As 
T — > T c from above, the scalar susceptibility averaged 
over '^-states' (rhs of (19)) diverges since a and n be- 
come degenerate. These modes are the analogue of the 
ones originally discussed by Hatsuda and Kunihiro Jjl[] 



using an effective model of QCD. For the present case, 
the calculation is done readily by rescaling $ — > y/V^ 
and noting that the quartic contribution in ( |l8| ) becomes 
subleading in large V. Hence, 



V \ \ 2 / 



1 



2go(T) 



(20) 



Near T c from above it is seen to diverge as 1/(T — T c ) 
with the critical exponent 7=1 (mean-field). We note 
that above the critical temperature, the level spacing is 
V^ 1 / 2 , which is intermediate between V^ 1 in the vacuum 
and y -1 / 4 in free space. Other sum rules are of course 
possible. 

The possibility of microscopic sum rules above the crit- 
ical temperature reflects on the assumption that the sym- 
metry restoration in two-flavor QCD is driven by uni- 
versality. The level correlations in the Dirac spectrum 
attests to the correlations still present on the 0(4) mani- 
fold despite the gap developing in the eigenvalue density. 
The gap is due to the fact that near zero virtuality the ac- 
cumulation of eigenvalues is not commensurate with the 
volume V. At high temperature, the quark eigenvalues 
are typically of order T, and both ir and a correlations 
are dissolved in the 'plasma' upsetting the present uni- 
versality arguments. In QCD we expect this to take place 
at a temperature T ~ 3T C [jl5| . 

6. The idea of microscopic sum rules above the crit- 
ical temperature may be readily checked using an (un- 
quenched) chiral random matrix model with one fla- 
vor. Again we may use the bosonized form of the parti- 
tion function (^J) , but now instead introduce the variable 
y = \Z~Nim ~ VNX, and rescale Q by Q — > Q = Qy/W. 
This leads to the following expression for the partition 
function 

poo 

Z{y) = t 2N / rdre^ 1 -^' J (2yr)e y2 = 
Jo 



Td 



t 2 



N t 2 - 1 



(21) 



Expanding the partition function in powers of y leads 
again to modified sum rules, the simplest example being 



1 



1 



\A fc >0 



NX 2 



t 2 -l 



(22) 



which is seen to diverge at t = 1 as expected. 

The comparison to numerical simulation is shown in 
Fig. H using a random set of matrices M distributed with 
a Gaussian measure. For high temperatures or large ma- 
trix sizes the agreement is good. At the critical tempera- 
ture, the finite size effects are important. Amusingly, we 
note the drop by two orders of magnitude at t ~ 3 = 3i c . 
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FIG. 2. The result ( |22[ ) (solid line) checked against numer- 
ical simulations using N = 20 (open squares), N = 50 (open 
circles) and N = 100 (open triangles) random matrices. 

The Green's function both for quenched (the mass play- 
ing the role of an external parameter) and one flavor chi- 
ral random matrix model in the rescaled variables may 
be readily obtained, 



G(y) 



1 



1 



y - 



NX, 



(23) 



The existence of a microscopic spectral density at 
T = T c for QCD opens up the interesting possibility 
of measuring both the critical temperature T c and the 
critical exponent 6 by simply monitoring the pertinently 
rescaled distribution z/»(s) of eigenvalues in lattice QCD 
simulation. We have also suggested that these correla- 
tions may persist even above T c . 
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Since for T > T c the eigenvalue spectrum develops a gap, 
there are no eigenvalues for small values of y, hence the 
resolvent is purely real. For finite sizes however, y may 
get out of the gap and our result breaks down as well 
as the scaling arguments. This is seen in Fig. ^ where 
for temperatures slightly above the critical one we are 
entering the nonzero eigenvalue density part of the spec- 
tra. This effect is shifted for higher values with increasing 
temperature and matrix size. 
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FIG. 3. Scaled and quenched resolvent (solid line) in 
comparison to a numerical simulation with N — 20 (long 
dashes), N = 50 (short dashes) and N = 100 (dotted line) 
random matrices. 



7. Using arguments based on universality we have sug- 
gested that the QCD Dirac spectrum may exhibit uni- 
versal spectral correlations at T = T c that reflect on the 
nature of the chirally restored phase. To illustrate our 
points, we have used a chiral random matrix model. Al- 
though the matrix model is based on a Gaussian weight, 
we provided physical arguments for why the results are 
insensitive to the choice of the weight at the critical point. 
In QCD this can be readily checked using our recent ar- 
guments j| at finite temperature, since the closest singu- 
larity to zero in the virtuality plane is persistently 'pionic' 
for T<T r [II. 
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